We have investigated the exclusive, radiative i?-meson decay to the 
charmless meson iC|(1430) in 88.5 x 10^ BB events accumulated between 
1999 and 2002 with the BABAR detector at the PEP-II storage ring. We 
present a preliminary measurement of the branching fractions B {B^ 
iC2*(1430)°7) = (1.22 ± 0.25 ± 0.11) x 10"^ and B {B+ K^{U30)+-f) = 
(1.44 lb 0.40 lb 0.13) X 10^^, where the first error is statistical and the second 
systematic. 
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Abstract 

We have investigated the exclusive, radiative S-meson decay to the cliarmless meson i^|(1430) 
in 88.5 X 10^ BB events accumulated between 1999 and 2002 with the BABAR detector at the 
PEP-II storage ring. We present a preliminary measurement of the branching fractions B{B^ — > 
K*{U30fj) = (1.22 ±0.25 ±0.11) x 10"^ and B{B^ /C*(1430)+7)= (1.44 ±0.40 ±0.13) x 10-^ 
where the first error is statistical and the second systematic. 
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1 INTRODUCTION 



In the Standard Model, flavor-changing neutral currents (FCNC) are forbidden at the tree level. 
For example, there is no direct coupling between the b quark and the s or d quarks. Effective 
FCNC are induced by loop (or "penguin") diagrams, where a quark emits and re-absorbs a W thus 
changing flavor twice, such as the b ^ transition depicted in Fig. ^ 




Figure 1: b —>■ sj "penguin" diagram. 

With large samples of i?-mesons and increasingly powerful background-suppression techniques, 
experimenters have succeeded in measuring Cabibbo-Kobayashi-Mashawa favored penguins. The- 
ory has also made good progress in computing the decay rates; in particular, recent completion 
of the NLO correction [l] to the inclusive 6 — > 57 decay rate has sparked experimental efforts to 
further improve the measurement. Comparisons between theoretical and experimental rates place 
strong constraints on physics beyond the Standard Model [2]. 

The discovery of -B ^ /C* (892)7 by the CLEO collaboration verified the existence of pen- 
guins. The same publication also reported the evidence for B (1430)7, later confirmed by 
the BELLE collaboration (4j. Measurements of the inclusive decay rate have been made by several 
experiments Detailed knowledge about the decays to resonant modes with masses higher than 
K*{892), such as ^ (1430)7 decay, will provide a better understanding of the inclusive 6 — > 57 
branching fraction in terms of the sum over exclusive modes. 

2 THE BABAR DETECTOR AND DATASET 

This study is based on 81.4 fb~^ of data collected at the T(45) resonance ( "on-resonance" ) with 
the BABAR detector at the PEP-II asymmetric e+ (3.1GeV)-e~ (9.0 GeV) storage ring, correspond- 
ing to 88.5 X 10^ BB pairs. We have also collected a data sample of 9.6 fb ^ at 40MeV below the 
T{4:S) energy ("off-resonance"). The number of BB meson pairs is determined from the ratio of 
the number of hadronic events to muon pairs in on- and off-resonance data. We assume that the 
T{4:S) decays equally to neutral and charged B-meson pairs. 

The BABAR detector is described elsewhere [H]. The analysis described below makes use of 
charged track, Ks and vr'' reconstruction, along with the charged particle identification. Charged 
particle trajectories are measured by a 5-layer double-sided silicon vertex tracker (SVT) and a 
40-layer drift chamber, which also provide ionization measurements (dE/dx) used for particle iden- 
tification. For charged tracks with momentum p > 1 GeV/ c, the measured transverse momentum 
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with respect to the beam axis (pt) has the resolution 



^ = 0.13%pT + 0.45%, (1) 

PT 

where pT is measured in GeV/ c. 

Photons and electrons are measured in the barrel and forward end-cap electromagnetic calorime- 
ter, consisting of 6580 Thallium-doped Csl crystals. The electromagnetic calorimeter resolution, 
aE, can be expressed as 

f = ^©1.9%, (2) 

where the energy E is measured in GeV. 

Charged particle identification is provided by the energy loss (dE/dx) in the tracking devices 
and by an internally reflecting ring-imaging Cherenkov detector (DIRC). The DIRC transports the 
Cherenkov light to a water- filled expansion volume equipped with approximately If, 000 photo- 
multiplier tubes. A K/tt separation better than four standard deviations is achieved for momenta 
below 3 GeV/c. 

We use Monte Carlo simulations of the BABAR detector based on GEANT 4 7_ to optimize 
our selection criteria and to determine signal efficiencies. These simulations take into account the 
varying detector conditions and beam backgrounds during the data-taking period. 



3 B CANDIDATE RECONSTRUCTION 

The /C2(1430) is reconstructed from K^, K^, and vr*^ candidates through the three modes 
Kf{U30) K+TT- and K^^iUSO) /sT+tt^, K°7r+. K° mesons are only reconstructed from the 
decay Kg vr+vr". In this paper the charge conjugate decays are implied unless otherwise stated. 

A photon candidate is defined as a localized energy maximum within the calorimeter acceptance 
—0.74 < cos^ < 0.93, where 9 is the polar angle to the detector axis. It must be isolated by 25 cm 
from any other neutral candidate or track and have a lateral energy profile consistent with a 
photon shower. To suppress photons from TT^{r]) decays, we veto any photon that combines with 
another photon of energy greater than 50 (250) MeV to form a 77 invariant mass in the range 
115(508) < M^^ < 155 (588) MeV/c^. 

The vr'^ candidates are reconstructed from pairs of photons, which have energy above 50 MeV 
and opening angle less than 36 degrees; the invariant mass of the two photons is required to be in 
the range 115 < Af^^ < 150 MeV/ c^; the vr*^ momentum is recalculated with a vr'^ mass constraint 
to improve the energy resolution. 

The and tt^ track candidates are required to be well reconstructed in the drift chamber and 
to originate from a vertex consistent with the e~^e~ interaction point (IP). A track is identified as a 
kaon if it is projected to pass through the fiducial volume of the DIRC, and the cone of Cherenkov 
light is consistent in time and angle with a kaon of the measured track momentum. A charged pion 
is identified as a track that does not satisfy the criteria for a kaon or an electron. 

The K'^ candidates are reconstructed from two oppositely charged tracks coming from a common 
vertex displaced from the IP by at least 0.2 cm. Each candidate is required to lie in a direction 
consistent with the momentum, while having an invariant mass within a 489 < M^+^~ < 
507 MeV/ window around the nominal K'^ mass. 
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The -ftr|(1430) candidate is required to have a Kir invariant mass within 120 and llOMeV/c^ of 
the known i^l*^ and K2~^ mass |0. For the K~^'ir~ mode, we require that the two tracks originate 
from a common vertex with a probabihty of at least 0.001. 

To remove Bhabha, radiative-Bhabha and tau events, we require that the second-order Fox- 
Wolfram moment [S] of the event is less than 0.9. 

The B candidates are reconstructed by combining one K2{14:30) and one 7 candidate with 
energy 1.80 < E* < 2.75 GeV in the rest frame of the T{4S). In this frame, the energy of the 
S-meson is given by the beam energy, ■v/s/2, which is known much more precisely than the energy 
of the B candidate. Therefore, to isolate the i?-meson signal, we use two kinematic variables: 
the difference between the reconstructed energy of the B candidate and the beam energy in the 
center-of-mass frame (A£^), and the beam energy substituted mass (m-Es), defined as 

^ES - V^boam -^'B' 

where 

-2'beam = a/s/2, 

Pb =Pk* +Pt 

For the modes containing a single photon candidate, namely K^n^ and K^n^, we adopt a 
technique from the CLEO analysis , which rescales the measured photon energy E* with a factor 
K, determined for each event, such that E'^, + kE* — E^^^^ = 0; this improves the tties resolution 
from 3.0 to 2.7MeY/c^. If we find multiple candidates with |A^| < 0.3 GeV and mss > 5.2 GeV/ 
in the same event, which happens in 3.1%, 6.3% and 4.9% of the events for the K~^7t~, -^5^^ 
and K~^TT^ modes, respectively, we take the candidate that has the invariant mass closest to the 
K2*(1430) mass jH|. 

The background has two components, one of which is continuum qq production, where q can be 
a u, d, s or c quark, with the high-energy photon originating from initial-state radiation or from vr*^ 
and r] decays. (This is known as the "non peaking" background in the m^s ^-^d distributions.) 

The second background contribution is from other B Xg^ modes (dominantly B — > if * (1410)7) 
and non resonant B Kirj decays. We label these the "peaking" background, since these decays 
have TTT-Es and AE distributions similar to the signal. 

In order to distinguish the B — > ii'g (1430)7 signal from B ii'* (1410)7 resonant 
decays, we examine the helicity angle distributions. The helicity angle 9h is defined as the angle of 
the or in the rest frame of the -fir|(1430) with respect to the flight direction of the (1^30). 
All three modes have different helicity angle distributions: syc^Ohcos^Oh for 7^2(1430), siv^Oh for 
if* (1410) and primarily siT?9H for non resonant decays assuming J = 1 for the spin of the Ktt 
system. The non resonant decays may have higher angular momentum contributions but the lowest 
possible angular momentum state is dominant, so the helicity angle distribution for non resonant 
decay is attributed to be the same as that oi B ^ if* (1410)7 decay. The systematic uncertainty 
related with this modeling is studied and included in the measured branching fraction uncertainty. 

4 BACKGROUND SUPPRESSION 

We have exploited the difference in the event topology between signal and continuum back- 
ground to reduce the continuum contribution, as well as the combinatorial BB background. The 
distribution of the thrust angle, defined as the angle between the direction of the photon candidate 
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and the thrust axis of the rest of the event (the particles not used in the reconstruction of the 
B candidate) in the center-of-mass frame, is shown in Fig. [21 We train a neural network with a 
combination of the thrust angle, the angle of the i?-meson candidate's direction with respect to the 
beam axis, the center-of-mass momentum vectors for tracks and photons not associated with the B 
candidate (binned in 10° intervals ranging from parallel to anti parallel relative to the photon mo- 
mentum), sphericity and the ratio of second-to-zeroth order Fox- Wolfram moments in the photon 
recoil system. The neural network improves background suppression significantly. The distribution 
of the neural network output is shown in Fig. |21 for MC signal, MC continuum background and 
off-resonance data. The neural network software used for this analysis is based on the Stuttgart 
Neural Network Simulator [lOj . 




Figure 2: Thrust angle distribution: Filled circles correspond to S — > -R'|(1430)7 signal MC 
simulation, the histogram corresponds to off-resonance data, and the open circles correspond to 
continuum background MC simulation. The vertical line indicates the cut value. 

The cuts on thrust angle (0t) and neural network output (NNO) have been optimized using 
an iterative method to minimize correlations. The final cuts are |cos0t| < 0.95 and NNO> 0.55, 
shown in Figs. |2l and El Because the variables used for the neural network training are mostly 
calculated from the rest-of-the-event information, we use a sample of fully reconstructed B — > D7r~ 
candidates in data, as well as a sample of simulated B Dtt~ events as control samples. The 
bachelor pion in the B Dtt~ decay is treated like the photon in — > i('|(1430)7 decay for the 
calculations of the event variables; the difference in the efficiency of the cut on the neural network 
output between different samples is taken as the systematic uncertainty related to this cut. 
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Figure 3: Neural network output distribution: Filled circles correspond to i? ^ K|(1430)7 signal 
MC simulation, the histogram corresponds to off-resonance data, and the open circles correspond 
to continuum background MC simulation. The vertical line indicates the cut value. 

5 ANALYSIS METHOD 

An unbinned maximum- likelihood technique is used to fit simultaneously the m^s, AE and 
I cos9h\ distributions. The fit is performed independently for each of the decay modes considered 
here. 

The signal ?ties and AE distributions are well described by an asymmetric resolution function 
("Crystal-Ball" function) 11_ S, having an approximately Gaussian core plus a long tail due to the 
energy leakage from the calorimeter for the photon candidates: 

{exp(— (m — mo)^/ (2(T^)) for m > mo — aa, 
(n/ay'exp(-a'^/2) r ^ (^) 

,, \! — / \n for m ^ ran — aa. 

The continuum background is parameterized empirically by a threshold function ^21 for m-ES and 
a linear function for AE. 

The cos 9h distribution of the signal has been parameterized with sin^^jycos^^jy — X{cos^9h — 
cos^0h), where A is a parameter determined from the Monte Carlo sample to account for the 
effect of the detection acceptance and efficiency. The | cos 9h \ distribution of the "non peaking" 
background is parameterized by the combination of exponential and constant components. 

Figures. 131 IS] and ini show the m-ESj and |cos0/f| distributions for the K^tt^, K^tt^ and 
K^TT^ modes, respectively. The \cos9h\ distributions in the signal region (defined as —0.15 < 
AE < 0.10 GeV, 5.272 < ttt-es < 5.288 GeV/c^) are also shown. The signal as well as background 
yields are allowed to vary in the fit. All the non peaking background parameters are allowed to 
float. The signal and peaking-background helicity angle, Crystal-Ball width and shape parameters 
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are constrained according to the MC expectations. The means of the signal rriEs functions 
are constrained to the MC expectations, cahbrated using B — > i^* (892)7 candidates from MC 
simulation and data, while the peaking-background means are allowed to float due to its complex 
composition. The signal and peaking-background yields are given in Tabled The signal significance 
has been evaluated from the change in the likelihood when the fit is repeated with the signal yield 
set to zero. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
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Figure 4: A fit to the (a) rriES, (b) A£; and (c) | cos 61/^1 distributions for B° i^|(1430)°7, 
ir|(1430)'^ — > K'^TT~ candidates in data, and (d) the | cos9h\ distribution in the signal region. The 
solid line shows a fit to the data. The peaking (dash-dot) and non peaking (dash) background 
contributions are also shown. 

Figure. 13 shows the Kit invariant mass distribution, which is fit with a relativistic Breit-Wigner 
function plus a first-order polynomial background, though we understand that it is not the most 
accurate description of the resonance due to possible interferences. There is a clear enhancement 
around 1.4GeV/c^ in the neutral mode and a slight enhancement in the charged modes. 

6 SYSTEMATIC STUDIES 

The total systematic error is the sum in quadrature of the components shown in Table |2j The AE 
resolution is dominated by the photon energy resolution, which is determined from data using vr^ 
and rj meson decays with symmetric decay-photon energies. The deviation in the reconstructed rj 
mass from the nominal rj mass provides an estimate of the uncertainty in the measured single photon 
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Table 1: The efficiency, fitted signal yield, peaking-background yield and measured branching 
fraction B{B K^{U30)-f) for each K|(1430) decay mode. 



Mode 


Efficiency 

(%) 


Signal 


Peaking 
background 


Signal 
significance 


B{B i^2* (1430)7) 
(10-5) 




6.4 


69. ± 14. 


47. ± 12. 


5.8 


1.22 ±0.25 ±0.11 




1.9 


29. ± 10. 


11. ±8. 


3.3 


1.69 ±0.59 ±0.20 




1.9 


20. ± 9. 


4. ± 7. 


2.6 


1.23 ±0.55 ±0.12 



energy. The photon isolation and tt^/t] veto efficiency depends on the event multiplicity and the 
effect is estimated by "embedding" MC-generated photons into both an exclusively reconstructed B- 
meson data sample and a generic S-meson MC sample. The photon and efficiency uncertainties 
are determined from a comparison of the efficiencies in data and MC for e~^e~ t^t~ events. The 
tracking-efficiency uncertainty is estimated from a sample of tracks well measured in the SVT. We 
estimate the uncertainties in the efficiency by comparing the data and MC distributions of the 
momentum and flight distance. The efficiency for kaon identification in the DIRC is derived from 
a sample of the decays D*~^ D^7r~^, with — > K~tt~^. 

The uncertainties in the background suppression cuts (the thrust angle and neural network cuts) 
are obtained by comparing the cut efficiencies between the B ^2 (1430)7 MC, and B — > Dvr"*" 
data and MC samples, and choosing the largest discrepancy. The invariant mass and width used 
for ^"2(1430) are from the PDG 0. The systematic error is obtained by varying the mean and 
width, within their errors, and using that which results in the largest change in efficiency. 

We estimate the systematic error due to the fitting procedures as follows. For the shape param- 
eters for m-ES) and | cos9h\ distributions, we vary the parameters in the fit within their errors 
from the MC expectations. We also test the validity of the peaking-background | cos Oh\ probability 
density function (PDF) by replacing it with different parameterized PDFs from MC samples. We 
double the largest deviation in these tests as the systematic error of the signal yield. There is also 
a small systematic error associated with the limited statistics of the signal MC sample. 

7 CONCLUSIONS 

We have presented a preliminary measurement of the branching fraction for B^ — > (1430)^7 of 
(1.22 ± 0.25 ± 0.11) X 10"^ with a 5.8cr statistical significance, which is an improvement over 
previous experimental results, while still being in agreement. We observe a signal with a statistical 
significance of 4.1cr for B~^ — > iir|(1430)''"7 and measure the preliminary branching fraction to be 
(1.44±0.40±0.13) X 10"^, by combining the results from Kg7r~^ and Ktt^ modes. Both results agree 
with the theoretical prediction from the relativistic form- factor model of Veseli and Olsson |13j . 
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Figure 5: A fit to the (a) toes, (b) and (c) \ cos6h\ distributions for S+ i^|(1430)+7, 
i^2^(1430) — > KgTr~^ candidates in data, and (d) the | cos distribution in the signal region. The 
solid line shows a fit to the data. The peaking (dash-dot) and non peaking (dash) background 
contributions are also shown. 
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Figure 7: K^ir (left) and K^tt^ and K^tt^ (right) invariant mass distributions for the signal 
region (see text for a definition of the signal region). 

Table 2: Fractional systematic uncertainties (%) in the measurement of B{B A'|(1430)7). 



Uncertainty 










S-counting 


1.1 


1.1 


1.1 


Photon detection efficiency 


1.3 


1.3 


1.3 


Photon energy scale 


1.0 


1.0 


1.0 


Photon energy resolution 


2.5 


2.5 


2.5 


Photon isolation 


2.0 


2.0 


2.0 


TT^/rj veto 


1.0 


1.0 


1.0 


/tt^ tracking 


1.6 


0.8 


0.8 


efficiency 




3.0 




TT^ efficiency 






2.5 


Sub-mode branching fraction 


2.4 


2.4 


2.4 


K|(1430) mass/width 


1.6 


1.0 


1.1 


Signal PDF parameters 


3.9 


5.8 


6.3 


Background suppression 


5.0 


6.0 


6.0 


Peaking-background modeling 


3.0 


5.0 


5.0 


MC statistics 


2.5 


3.2 


3.2 


Total 


9. 


12. 


10. 



[12] The background for each mode is fit with dN/dniEs oc ■Mes_ ^ , /i _ ^^es xexp(— ^(1— -y/-^ )), 

■C-beam y ^beam ^beam 
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